Damage and Healing in Fatigue Fracture 
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We present an experimental and theoretical study of the fatigue failure of asphalt under cyclic 
compression. Varying the load amplitude, experiments reveal a finite fatigue limit below which 
the specimen does not break, while approaching the tensile strength of the material a rapid failure 
occurs. In the intermediate load range, the lifetime decreases with the load as a power law. We 
introduce two novel theoretical approaches, namely, a fiber bundle model and a fuse model, and show 
that both capture the major microscopic mechanisms of the fatigue failure of asphalt, providing an 
excellent agreement with the experimental findings. Model calculations show that the competition 
of damage accumulation and healing of microcracks gives rise to novel scaling laws for fatigue failure. 

PACS numbers: 62.20.Mk, 46.50.+a, 61.82.Pv 



The fracture of disordered media represents an impor- 
tant applied problem, with intriguing theoretical aspects. 
Statistical models have been successfully applied in the 
past to analyze fracture under quasistatic conditions, but 
the effect of cyclic loading is less explored . Laboratory 
experiments reveal that fatigue failure under repeated 
loading is due to a combination of several mechanisms, 
among which damage growth, relaxation due to viscoelas- 
ticity, and healing of microcracks play an essential role 
CJ 13 Theoretical approaches have serious difficul- 
ties to capture all these mechanisms 0, IE S S IS 
and fatigue life prediction is still very much an empiri- 
cal science. Understanding this problem has crucial im- 
plications even for everyday applications. For example, 
fatigue failure occurring in roads due to repeated traf- 





FIG. 1: Setup of the experiments, a) A cylindrical asphalt 
sample is subjected to diametrical compression applied pe- 
riodically. FBM discretizes the region where tensile stress 
emerges (white rectangle) in terms of fibers, b) At complete 
failure a crack spans the cylinder along the load direction. 
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fic loading cause main distress, limiting the lifetime of 
asphalt pavements. 

In this Letter we present a detailed experimental and 
theoretical study of the fatigue performance of hot mix 
asphalt (HMA). We carried out fatigue life tests of spec- 
imens measuring the accumulation of deformation with 
the number of loading cycles and the lifetime of speci- 
mens varying the load amplitude. To obtain a theoretical 
understanding of the experimental findings, we worked 
out two novel modelling approaches for fatigue failure, 
namely, a fiber bundle model 0, and a fuse model 
|Toj| . We then show that both descriptions capture the 
stochastic nature of the fracture process, the immediate 
breaking of material elements and the cumulative effect 
of the loading history. Two physical mechanisms are con- 
sidered which limit the accumulation of damage: a finite 
activation threshold of crack nuclcation below which the 
local load does not contribute to the ageing of the mate- 
rial and healing of microcracks under compression, which 
leads to damage recovery. The analytical and numerical 
results of the model calculations provide a good quan- 
titative agreement with the experimental findings. We 
show that the competition of nucleation and healing of 
microcracks leads to novel type of scaling laws for fatigue 
fracture with universal scaling exponents. 

In order to obtain a quantitative characterization of the 
process of fatigue failure, we carried out fatigue life tests 
of asphalt under cyclic diametric compression of cylindri- 
cal specimens at a constant external load ctq (sec Fig. ^| . 
HMA is the primary material used to construct and main- 
tain pavements and roadways due to its good mechani- 
cal performance and high durability. From the structural 
point of view asphalt is a combination of aggregates (usu- 
ally crushed stone and sand), filler (cement, hydrated 
lime or stone dust) and a bituminous binder. Cylindri- 
cal samples of HMA were produced using the Marshall 
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FIG. 2: Deformation as a function of the number of loading 
cycles. The continuous lines result from our theory. Inset: 
Load on the fibers in FBM as function of time at different val- 
ues of ao I o~ c for uniformly distributed threshold values setting 
t = oo in Eq. . 
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FIG. 3: ( Color online) Nf as function of load ao/ac for FBM 
varying the value of r. FBM provides an excellent fit of the 
experimental results with 7 = 2.0, r = 15000, and a — 0.01. 
Fuse model results are presented in the inset varying the value 
of 7 and the threshold current io . 



method |ll| and then loaded by a hydraulic device. Un- 
der repeated loading at a constant amplitude ao, the de- 
formation e was monitored as a function of the number 
of cycles N cyc i e . Furthermore, the total number of cycles 
to complete failure Nf was measured varying oq. Figure 
El presents representative examples of s(N cyc i e ) recorded 
at loads 30% and 40% of the tensile strength a c of the 
specimen. It can be observed that due to the gradual ac- 
cumulation of damage, the deformation e caused by the 
same load <7q monotonically increases until catastrophic 
failure occurs after a finite number of cycles Nf. The 
derivative of e(N cyc i e ) also shows a monotonous increase 
and diverges when approaching the point of macroscopic 
failure. Increasing the external load the functional form 
of e(N cyc i e ) remains the same, however, the lifetime of 
the specimen Nf gets shorter. The fatigue lifetime Nf 
measured at different fractions of the tensile strength a c 
(Fig. 13 reveals the existence of three distinct regimes. 
First, approaching the tensile strength of the material 
< J o/< J c —> 1 the lifetime Nf rapidly decreases indicating 
an immediate failure of the specimen. At the other ex- 
treme, a lower threshold value of the external load 
can be identified below which the specimen suffers only 
partial damage giving rise to an infinite lifetime (fatigue 
limit). In the intermediate regime the experimental re- 
sults follow a power law known as Basquin law 0,00,0] 

»,~(?y. (i) 

where a = 2.2 ± 0.1, as shown in Fig.|3 

The experiments show that the fatigue crack growth is 
localized to a narrow region between the loading plates 



(see Fig. |2>) where locally a tensile stress emerges per- 
pendicular to the external load . To give a theoretical 
description of the failure process, we focus on this region 
and discretize it by a fiber bundle model (FBM) as illus- 
trated in Fig. H 0] ■ We consider a bundle of parallel 
linear elastic fibers with the same Young modulus E. Un- 
der diametrical compression of the disc-shaped specimen, 
the fibers experience a tensile loading and gradually fail 
due to immediate breaking or to the ageing of material 
elements 0. More precisely, the following two mecha- 
nisms are considered: (I) fiber i (i = 1,...,N) breaks 
instantaneously at time t when its local load Pi(t) ex- 
ceeds the tensile strength p\ h of the fiber. (II) All intact 
fibers undergo a damage accumulation process due to the 
load they have experienced. The amount of damage Ac$ 
occurred under the load Pi(t) in a time interval At is 
assumed to have the form Aci = api(t) J At hence, 
the total accumulated damage Cj(i) until time t can be 
obtained by integrating over the entire loading history 
[lflj . The exponent 7 > controls the damage accumu- 
lation rate and a > is a scale parameter. The fibers 
can only tolerate a finite amount of damage and break 
when Ci(t) exceeds a threshold value c\ h . The two break- 
ing thresholds p l th and c\ h are random variables with a 
joint probability density function h(j>th,Cth)- Assuming 
independence of the two breaking modes, h can be fac- 
torized into a product h(p th ,c t h) = f(cth)g(pth), where 
f(cth) and g(pth) are the probability densities and F(cth) 
and G(pth) the cumulative distributions of the breaking 
thresholds pth and Cth, respectively. For simplicity, af- 
ter each breaking event the load of the broken fiber is 
equally redistributed over the intact ones irrespective of 
their distance from the failure point (global load sharing) 



3 






□ 


unifoi'm, 1= 1.0 ^ 






unifoi'jn, 7= 2J3 i 




WeibuU, j= 1.0 ' 






Weibull, t= 0.5 ; 









tr to c 



FIG. 4: (Color online) a) Deformation rate dp/dt as func- 
tion of the time to macroscopic failure for different disorder 
distributions and 7 = 1. Straight lines are drawn to guide 
the eye with slopes 1.5 and 1.8. b) tf as a function of o~o/a c 
for uniform and Weibull distributions (m = 2.0) varying the 
value of 7. The slope of the straight lines is equal to the value 
of the corresponding exponent 7. 



Under a constant tensile load <tq, the load on a single 
fiber po is initially determined by the quasi-static consti- 
tutive equation of FBM a = [1 - G(p ))p H 0. The 
external load ctq must fall below the tensile strength of the 
bundle (To < cr c , otherwise the bundle will fail immedi- 
ately. As time elapses, the fibers accumulate damage and 
break due to their finite damage tolerance. These break- 
ings, however, increase the load on the remaining intact 
fibers which in turn induce again immediate breakings. 
This way, in spite of the independence of the threshold 
values pth and c t h, the two breaking modes are dynam- 
ically coupled, gradually driving the system to macro- 
scopic failure in a finite time tf at any load values (Tq. 
Healing of microcracks can be captured in the model by 
introducing a finite range r for the memory, over which 
the loading history contributes to the accumulated dam- 
age mm. 

Finally, the evolution equation of the system 
can be cast in the form 



a = [1 - F(a / e 



(*-*') 



p{t'Vdt')] [1-G(p(t))]p(t), 



(2) 

where the integral in the argument of F provides the 
accumulated damaged at time t taking into account the 
finite range of memory by the exponential term [Hj . In 
principle, the range of memory r can take any positive 
value r > such that during the time evolution of the 
bundle the damage accumulated during the time interval 
t' < t — t heals. Equation J2J) is an integral equation 
which has to be solved for the load p(t) on the intact 
fibers at a given external load <7o with the initial condition 
p(t = 0) = po- The product in Eq. © arises due to 
the independence of the two breaking thresholds. Wc 
note that Eq. 10) recovers the usual constitutive behavior 
of FBM 0, [i| when damage accumulation is suppressed 
cither by increasing the exponent 7 or decreasing the 



range of memory r — > 0. 

The inset of Fig. presents examples of the solution 
p(t) of Eq. (J2J obtained for breaking thresholds uniformly 
distributed in the interval [0, 1] at different ratios oo/cc 
setting t — > 00. Since p(t) is simply related to the 
macroscopic deformation e of the bundle p(t) = Ee(t), 
these results can directly be compared to the experimen- 
tal findings. The agreement seen in Fig. [21 is obtained 
using Weibull distributions P(x) = 1 — exp [— (x/X) m ] 
for the two breaking thresholds in Eq. J5J with the same 
Weibull exponent m and scale parameter A. Our calcu- 
lations also reveal that approaching macroscopic failure, 
the derivative of p(t) has a power law divergence as a 
function of the time to failure h? ~ (tf — t)~P . Fig. 0Ji 
shows that the exponent j3 solely depends on the type of 
disorder. Specifically, for breaking thresholds distributed 
over a finite and infinite range, we obtain the exponents 
(3 = 1.5 ± 0.02 and f3 = 1.8 ± 0.06, respectively, defining 
two different universality classes of the fatigue failure. 

It is possible to recover the Basquin law Eq. (JJJ from 
Eq. Q, i.e. it can be shown analytically that for 
co/c c << 1 and r — > 00 the lifetime of the system has a 

power law dependence on the external load tf ~ (f 2 ') > 
where 7 is the damage accumulation exponent, indepen- 
dent on the type of disorder. Figure^ shows that the nu- 
merical results are in excellent agreement with the above 
analytic prediction. 

Due to Eq. (J2J, without healing (t — > 00) the cumu- 
lative effect of the loading history gives rise to a macro- 
scopic failure of the system at any load. However, our 
experiments revealed that damage recovery caused by 
healing of microcracks results in a finite fatigue limit 07 , 
below which the sample does not break. Since healing 
takes place in the polymer binder, it can be controlled by 
changing the temperature Q . In our FBM the healing of 
microcracks is captured by the finite range of memory r. 
This is illustrated in Fig. presenting the lifetime tf for 
a fixed load varying the value of r over a broad range. 
It can be seen that by decreasing r, the lifetime of the 
system increases and, due to the competition between nu- 
cleation of new microcracks and healing of the existing 
ones, a finite critical value r c emerges below which the 
system only suffers a partial damage and has an infinite 
lifetime tf — > 00. Plotting the lifetime tf as function of 
the distance from the critical point r c (inset of Fig.[SJ) we 
find a power law dependence of t f 



(3) 



where 5 = 0.5 ± 0.01 is an universal exponent, indepen- 
dent of the type of disorder and of the damage accumu- 
lation exponent 7. Consequently, at a given temperature 
where the system is characterized by a fixed value of r, 
a finite fatigue limit 07 emerges at which (Fig. the 
lifetime diverges 

Besides healing, another important mechanism which 
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FIG. 5: Lifetime tf of FBM at a constant load ao/a c = 0.7 
varying r. Approaching r c the lifetime diverges. Inset: tf as 
a function of the distance from the critical range of memory 
r — r c for different disorder distributions and 7 exponents. 
The horizontal lines indicate the value of tf(r = 00). 



limits damage accumulation is a finite activation thresh- 
old of microcrack nucleation. In order to study this effect 
we consider the random fuse model (RFM) [lOj of frac- 
ture and extend it by introducing a history dependent 
ageing variable of fuses. We construct an L x L tilted 
square lattice of initially fully intact bonds with identi- 
cal conductance but random failure thresholds i c . The 
threshold values i c are uniformly distributed between a 
small current value io and 1 (io << 1). For a given value 
of current / applied between two bus bars of the lat- 
tice, the local current through each bond is determined 
by solving numerically the Kirchhoff equations. Fuses 
burn out irreversibly when the current exceeds the lo- 
cal failure thresholds. This process is then followed by 
the recalculation of the current values. In order to cap- 
ture fatigue cracking in the model, intact fuses are as- 
sumed to undergo an ageing process, modeled by a vari- 
able A{t) = 53f=i a (*(^') — *o) 7 - A fuse fails due to fa- 
tigue when A(t) > A max , where A max is a failure thresh- 
old uniformly distributed between 1 — 6 and 1 + 6 with 
6 = 0.1. Comparing the two modelling approaches, the 
ageing variable A(t) of RFM is analogous to the accumu- 
lated damage c(t) of FBM, however, only current values 
above io contribute to A(t), which captures the finite ac- 
tivation threshold of microcrack nucleation. The inset of 
Fig. 13 demonstrates that RFM of ageing fuses provides 
qualitatively the same behavior as FBM, i.e. rapid fail- 
ure at high current values /, a Basquin regime Eq. Q at 
intermediate currents with an exponent equal to 7 and a 
finite fatigue limit 07 determined by the threshold current 

Jo- 
in summary, our experiments on the fatigue failure of 
asphalt under cyclic compression revealed three regimes 



of the failure process depending on the load amplitude: 
instantaneous breaking, a Basquin regime of a power 
law decrease of lifetime and the existence of a fatigue 
limit below which no failure occurs. We introduced two 
novel modeling approaches both capturing the essential 
ingredients of the fatigue failure of bituminous materi- 
als. These models provide a comprehensive description of 
the experimental findings and additionally revealed novel 
scaling laws of fatigue fracture: approaching the macro- 
scopic failure, the process of fatigue fracture accelerates 
and is characterized by a finite time power law singular- 
ity of the deformation rate. The exponent is different for 
bounded and unbounded disorder distributions defining 
two universality classes of fatigue fracture. Due to the 
interplay between damage and healing, at each load level 
a critical value of the range of memory emerges where 
the lifetime of the system has a universal power law di- 
vergence. 
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